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Abstract 

t— s For the problem of diffraction of harmonic scalar waves by a lossless periodic slab 

scatterer, we analyze field sensitivity with respect to the material coefficients of the 
slab. The governing equation is the Helmholtz equation, which describes acoustic or 

i— i electromagnetic fields. The main theorem establishes the variational (Frechet) deriva- 

five of the scattered field measured in the H 1 (root-mean-square-gradient) norm as a 
function of the material coefficients measured in an L p (p-power integral) norm, with 
2 < p < do, as long as these coefficients are bounded above and below by positive 
constants and do not admit resonance. The derivative is Lipschitz continuous. We 
also establish the variational derivative of the transmitted energy with respect to the 
material coefficients in L p . 

CO 

Key words: Periodic slab; open waveguide; guided modes; scattering; Helmholtz equation; 
variational calculus; transmission coefficient; sensitivity analysis; elliptic regularity. 

OO 
CN 

cN 1 Introduction 



OO 

o 



This work treats the variational calculus of time-harmonic fields scattered by a periodic slab 
structure as functions of the material coefficients of the scatterer. We deal with scalar fields u 
governed by the linear Helmholtz equation 



X 

S V-tVu + cAm = 



which governs acoustic fields and, in case the coefficients are invariant in one direction, po- 
larized electromagnetic fields, in a composite material characterized by the spatially varying 
coefficients e and r. We take these coefficients to be real and positive, which means that 
the structure is lossless. Figure [T] depicts an example of the type of scatterer we consider. 
The slab is periodic in two directions and finite in the other, and it is in contact with the 
ambient space, making it an open waveguide. A traveling time-harmonic wave, originating 
from sources exterior to the slab, is incident upon the slab at an angle and is diffracted by it. 
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Figure 1: An example of a periodic slab scatterer. The slab is infinite and periodic in 
the x\ and X2 directions (only sixteen periods are shown) but of finite thickness in the X3 
direction. We impose no fixed boundary conditions (as Dirichlet or Neumann); rather, the 
slab is in natural contact with the ambient space to its left and right, which is homogeneous 
with e(x) = £0 and t(x) = To- This figure depicts a slab consisting of two homogeneous 
materials. In this work, we treat the much more general case, in which e{x) and t(x) are 
measurable functions bounded from below and above. The reflected and transmitted fields 
depict the diffractive orders associated with the angle of incidence of the source field. 



Our aim is to compute the sensitivity of the resulting total field to variations of the material 
properties (e and r) and geometry of the slab, as well as the sensitivity of the amount of 
energy transmitted across the slab. 

A motivation for this subject is the desire to optimize the way in which energy flows 
through a periodic slab or film, as well as the related inverse problem, in which one seeks 
to determine the structure that produces given diffracted field patterns upon illumination 
by plane waves. Slabs of photonic crystal structures can be used to guide energy of an 
incident wave at specific frequencies through channels to the other side of the slab [18J. 
The characteristics that one seeks to optimize are the amount of transmitted energy and the 
directionality of the field that is transmitted, as well the electromagnetic mode density, which 
is important for control of the spontaneous emission rate of atoms placed in the structure [3] • 
The variational calculus of the scattered, or diffracted, field as a function of the structural 
parameters is the basis for control and optimization of these properties. 

Variations of practical interest are not typically uniformly small across the scatterer; 
rather they tend to be large but supported in a small domain. For example, one may wish 
to vary the diameter of a dielectric sphere S of fixed permittivity S\ repeated in a two- 
dimensional periodic array within a matrix of permittivity Eq or the diameters of the holes 
in the example of Figure [TJ The function e = £\Xs + £ oXs c (Xs * s the characteristic function 
of S) is not continuous with respect to the diameter of S if the function is measured in 
the supremum norm, or L°° norm, but it is continuous if the function is measured in any 
p-power integral norm, or L p norm, with 1 < p < 00. Therefore, the question of whether the 
scattered field is differentiable, or even merely continuous, with respect to IP perturbations 



of the scatterer is an important one. In this work, a rigorous formulation (Theorems 15 
and 16) of the following theorem is proved: 
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Theorem The scattered field of a lossless periodic slab as well as the transmitted energy, 
for a fixed incident wave, are Frechet differentiable with respect to the coefficients e and t, 
if the field and its gradient are measured in the root-mean-square norm (Sobolev norm H 1 ) 
and e and t are measured in an LP norm, with p < oo, as long as e and r are bounded from 
above and below and do not admit resonance. The derivatives are Lipschitz continuous. 

Frechet differentiability with respect to the material coefficients in an LP norm implies 
Holder continuity with respect to variations of a smooth boundary of a homogeneous com- 
ponent of the scatterer. In fact, the scattered field has been shown to be differentiable with 
respect to variation of periodic interfaces separating materials of differing dielectric coeffi- 
cient in two-dimensional polarized electromagnetic scattering problems. See, for example, 
Bao [2], Dobson [7], and Elschner and Schmidt (TTJ [TU], as well as [El E], [Z], and Bao and 
Bonnetier pQ for applications to optimal design. The differentiability of solutions to strongly 
elliptic equations in a bounded domain as well as functionals of these solutions, with respect 
to the boundary in norms of Holder continuity, is treated by Pironneau [21] (§1.7, 6.2). In 
their study of the inverse problem for bounded impenetrable obstacles, Colton and Kress [5] 
(§5.3) prove the Frechet differentiability of the far field pattern in the the L 2 norm of the 
sphere as a function of the boundary in the norm of continuous differentiability. The inverse 
problem for scattering by periodic interfaces is treated by Kirsch [13] and in [TT] . 

Theorem [0] implies differentiability with respect to any L q norm with q > p. Obtaining an 
upper bound on the minimal p is an open problem, whose solution would facilitate numerical 
implementation of the variational gradient. The formal calculus of variations leads to a 
candidate for the gradient of the field and transmission coefficient as functions of e and r. 
The gradient of the transmission coefficient is expressed in terms of an adjoint problem, 
derived formally by Lipton, Shipman, and Venakides [16]. In that work, the authors used the 
formal results in a two-dimensional reduction (where e and r are constant in one direction) 
to manipulate numerically the transmission coefficient as a function of frequency by varying 
the slab structure. In this work, this gradient is established rigorously for L p perturbations 
of e and r. 

The proof uses N. Meyers' theorem on higher integral regularity (p > 2) of solutions 
of elliptic equations and their gradients |17j . In order to apply the theorem, one needs an 
a priori bound on the solution of the scattering problem that is independent of the material 
coefficients. The obstruction to such a bound is field resonance in the structure, resulting 
from the presence of guided modes. A guided mode is a pseudoperiodic solution (Bloch 
solution) to the Helmholtz equation that falls off exponentially with distance from the slab. 
Mathematically, it is self-sustained, that is, not forced by an incident source field. Because 
a solution of the scattering problem is not unique for a given structure at a frequency and 
Bloch wavevector that admit a guided mode, the scattered field is not uniformly bounded 
near these parameters. This work concerns the perturbation of the material properties within 
a range that excludes resonance, in which the scattering problem necessarily has a unique 
solution. Perturbation analysis near resonance is singular, and quite a different problem, as 
the field and transmitted energy exhibit anomalous behavior near a guided mode frequency. 
Rigorous perturbation analysis with respect to frequency and Bloch wavevector about a 
guided mode is presented in [23], and similar analysis with respect to material coefficients 
and geometry is possible. 
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The exposition of the ideas and results is summarized as follows. 



Section[2]presents the mathematical formulation of the problem of scattering of incident 
traveling waves by a periodic slab structure as well as formal perturbation analysis. 
This analysis gives rise to the correct candidate for the derivative of the scattered 
field u with respect to variations of the material coefficients e and r. 

The sensitivity of the transmitted energy to variations of the scatterer is discussed in 
Section |3j with specialization to structures with homogeneous components. 

Section [4] develops the weak formulation of the scattering problem, in which the fre- 
quency and Bloch wavevector are parameters. We discuss eigenvalues of the sesquilin- 
ear forms associated with the scattering problem and their relation to guided modes 
of the slab. 

The main contributions of this work are stated and proved in Section [5] The first, 



Theorem 13 (p. 18), establishes an a priori bound on the root-mean-square norm of 
the solution of the scattering problem and its gradient in the scatterer as long as the 
structure, frequency, and wavevector do not admit a guided mode. This result, together 



with Meyers' regularity theorem are used to prove the main result, Theorem 15 (p. p3 



on the differentiability of the scattered field with respect to L p variations of the material 



coefficients. Theorem 16 applies the main theorem and the adjoint method to give 
an explicit representation of the variational gradient of the transmitted energy as a 
function of the material coefficients. 



2 The scattering problem and sensitivity analysis 

The aim of this section is to derive a candidate for the variational gradient of the field 
scattered by a periodic slab as a function of the material coefficients of the scatterer. The 
variational calculus is treated rigorously in Section [5j First, we present the mathematical 
formulation of the scattering problem. 



2.1 The scattering problem 

We shall consider time-harmonic solutions U(x,t) = u(x)e~ lut (x = (xi,X2,x^) G 1R 3 ) of the 
scalar wave equation 

^ = V-rVf/, (1) 

in which the material coefficients e(x) and r(x) are positive, 27r-periodic in x\ and X2, and 
bounded from below and above. The spatial factor u satisfies the Helmholtz equation 

V'tVm + £w 2 m = 0. (2) 

By means of the (partial) Floquet transform in (xi,x 2 ), a solution u can be decomposed 
into an integral superposition of components u(x;k), where k = (ki,/^) £ ^- 2 , that are 
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n-pseudoperiodic in x\ and X2 with periods 27r (see [H] or [23], for example). This means 
that u(x; k) satisfies 

V-tVm + elu 2 u = 0, (3) 
u(x; k) = u peI (x; K ) e 1 ( Klx i+ K2:c 2) anc j u ^ has period 27r in x\ and x 2 . (4) 

The Block wavevector k = k^) is related to the angle of incidence of an incoming wave 

e «r7a;3 e j((mi+fti)xi + (?Ti2+K2)3;2) ^ 

for some m = (m 1 ,m 2 ) G Z 2 , which impinges upon the left side of the slab at an angle of 
9 = arctan \m + k\/t] with the normal. We shall take k to lie in the first Brillouin zone, 

K ^ [ _ 2' I) 2 ' 



1 1^2 

2' 2' 



because each K6l 2 can be written k = m + R for m G Z 2 and k£ [- 

Exterior to the slab (£3 < z~ and £3 > z + ), where the material is homogeneous, we set 
e = £0 > and r = t > 0, and the periodic factor be decomposed into Fourier 

components parallel to the slab (in x' = (a^i,^)). They are indexed by m G Z 2 (with 
different coefficients on the two sides of the slab), and the x 3 -dependence of each component 
is determined by separation of variables in the equation V • Vu + u 2 (£ /t )u = 0, 

u(x',x 3 ; k) = (c£0+(s s ) + c-0-(x 3 )) (6) 

m£Z 2 

in which 0^ are independent solutions of the ordinary differential equation 0^ + rj^m = 0, 
where the numbers r\ m are defined through 

1]^ + \m + k\ 2 = uj 2 £ /t . (7) 

The Fourier harmonics are known as the diffractive orders or diffraction orders associated 
with the periodic structure. There are many references that expound these ideas, including 
Wilcox and Neviere [20] ■ The </>^j are either oscillatory, linear, or exponential, depending 
on the numbers r) m . We make the following definitions: 

m G Z p •<==>■ 77^ > 0, r\ m > (propagating harmonics), 

m G Zi -<=>- rj 2 m = 0, i] m = (linear harmonics), (8) 

m G Z e •<=>- 77^ < 0, — irjm > (evanescent harmonics). 

In the problem of scattering of source fields given by traveling waves impinging upon the 
slab, we must exclude exponential or linear growth of u as \x 3 \ — > 00. The form the of the 
total field is therefore (we suppress the K-dependence of u(x', x 3 ] k)) 

u(x',x 3 ) = aZe iVmX3 e i{m+K)x ' + ^ a m e~ iVmXs ■■ e ^ m+K)x ' {x 3 < (9) 
u(x',x 3 ) = bZe- lVmX3 e i{m+K)x ' + b m e iVmX3 e i{m+K)x ' {x 3 > z + ). (10) 
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The infinite series are understood in the L 2 sense. The first sums in these expressions 
represent right-traveling source waves incident upon the slab from the left side and left- 
traveling source waves incident upon the slab from the right side. We say that a function u 
is outgoing if it is of the form (9p0) with a™ c = and 6™ c = for all m G Z p . 



Definition 1 (Outgoing and incoming) A complex-valued function u defined on IR 3 is 
said to be outgoing if there are real numbers z_ and z + and sequences {a m }^oo an d {^m}-oo 
in £ 2 (— oo, oo) such that 



u[x) 

meZ 2 



J2 a m e- ir > mX3 e i( - m+K)x ' for x 3 < z., (11) 

nez 2 

u(x) = b m e irimX3 e i{m+K)x ' for x 3 > z+. (12) 

The function u is said to be incoming if it admits the expansions 

u{x) = a m e tVmX3 e t{m+K)x ' for x 3 < z_, (13) 



u(x) = 



b m e- iVmX3 e i{m+K)x ' forx 3 >z + . (14) 

We shall take the pseudoperiodic source field to be a superposition of traveling waves 
incident upon the slab from the left and right. We think of these waves as emanating from 
x 3 = — oo and from x 3 = oo: 

u hlc (x',x 3 ) = (< c e ir > mX3 + ^ e -iVmX 3 y{m+K)x' _ ( 15 ) 

The problem of scattering of the incident wave w inc by the slab is expressed as a system 
characterizing the total field u, which is the sum of the incident field u mc and the scattered, 
or diffracted, field u sc , the latter of which is outgoing. The "strong form" of the problem is 
posed for functions e and r that are smooth except on a set £ consisting of continuously 
differentiable surfaces of discontinuity, with normal vector n. The "weak form", presented 



in subsection 4.1, allows e and r to be merely measurable. 



Problem 2 (Scattering of an incident wave, strong form) Given an incident field (15), 
find a function u that satisfies the following conditions. 

V-rVu + uo 2 eu = znM 3 \S, (16) 

u and rd n u are continuous on E, (17) 

u is k -pseudoperiodic in (xi,X2), (18) 

u = u mc + u sc , with u sc outgoing. (19) 

One can generalize the scattering problem by introducing sources originating from the 
slab itself or from points outside the slab. Such sources are represented by a periodic function 
h and a periodic vector field £, which enter the equation thus: 

V- rVu + lu 2 eu = V-^ + h. (20) 



6 



In our investigation of the perturbation of the scattering Problem |2j we will be concerned 
with an auxiliary problem involving sources that are confined to the region between x 3 = z_ 
and x 3 = z + (besides the incident source field originating from x 3 = ±00). 



Problem 3 (General scattering, strong form) Given an incident field (15), find a func- 
tion u that satisfies the following conditions. 

V-rVw + co 2 £u = V-£ + h for z^ < x 3 < z + and x S, (21) 

V-rV« + w 2 e« = otherwise, (22) 

u and rd n u are continuous on S ; (23) 

u is n-pseudoperiodic in (xx,x 2 ), (24) 

u = u mc + u sc , with u sc outgoing. (25) 

2.2 Formal sensitivity analysis 

Let u be the solution of the scattering Problem [2] (existence and uniqueness will be dealt 
with later), and let u + u be the solution of the scattering problem with the same incident 
field but with e + e and r + f in place of e and r. The coefficients e and r exterior to the 
slab remain fixed. The functions u and u + u satisfy 

V-tVu + uj 2 eu = 0, (26) 
V- (r + f)V(u + u) + uj 2 (e + e)(u + u) = 0. (27) 

Subtracting these equations yields the equation for the perturbed field it, 

V- (r + r)V« + u 2 (e + e)« = -V- rVn + u 2 eu, (28) 

and u is outgoing because the incident fields for u and u + u are identical and the forcing 



term on the right-hand side of (28) is confined to z_ < X3 < z + . If we remove the terms on 
the left-hand side that are quadratic in e, f, and u, we obtain the differential equation for 
the formal leading-order sensitivity u of the total field as a function of the perturbations 
e and r. We denote this linear approximation to u by u , 

V- tVu + u 2 eu = -V- fVu + u 2 eu, (29) 
Mo is outgoing. (30) 

In order to establish that the linear map (e, f) 1— > uq is truly the variational differential 
of u with respect to (e, r), we should demonstrate two things, 

\\u \\<C\\(s,t)\\, (31) 
11 """ ^Oas ||(e,f)|| ->0. (32) 



The appropriate norm in which to measure w is the i7 1 -norm, restricted to a domain Q 
(Figure 2) comprising one period of the structure between X3 = z_ and X3 = z + : 

Q = {x G M 3 : < xi < 2tt, < x 2 < 2n, < x 3 < z + }, (33) 

r± = {x e M 3 : < xx < 2tt, < x 2 < In, x 3 = z±}. (34) 
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Figure 2: The rectangular prism $7 encloses one period of the slab structure, which is defined 
by material coefficients e(x) and r(x) for x = (x\, X2, x 3 ) G f2 that are extended periodically 
in the xi and X2 directions. For x 3 < z_ and for X3 > z+, the material is homogeneous, 
with e(x) = Eq and t(x) = tq. 



The two-dimensional squares T± are the left and right boundaries of f2. The normal vector 
n to T is taken to be directed outward, so that 



d n u 



—du/dx 3 on T_, 
du/dx 3 on r + . 



(35) 



The H 1 norm in Q is 



r \ 1/2 

J (\Vu\ 2 + \u\ 2 )dV\ . 



(36) 



The main result of this work proves that (31 ) and (32) hold if (e, r) is measured in some LP 
norm in Q, with p < 00, 



{£,T)\\ L p(n) 



[ (\e\ p + \f\ p )dV 
in 



i/p 



(37) 



3 Energy transmission and special structures 

We apply the results discussed in the previous section to the sensitivity analysis of the 
amount of energy of an incident wave that is transmitted from one side of the slab to the 
other. The formal analysis of the adjoint problem associated with the differential of the 
transmitted energy that was derived in [16] is revisited in the light of the rigorous results of 
this work. 



3.1 Variation of the transmitted energy 

Let us send a traveling wave toward the slab from the left and consider the energy transmitted 
to right side of the slab. This means that we take 6™ c = in (10). We are interested in the 
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sensitivity of the transmitted energy to perturbations of the material coefficients e and r. 
The time-averaged energy flux through one period of the right-hand boundary of the slab is 
defined by 

£ = Im / r ud n u, (38) 



in which u is the solution to the scattering Problem [2] This quantity can be expressed in 
terms of the Fourier coefficients b m of the propagating harmonics of the transmitted field: 

u = b m e iVmX3 e i{m+K)x ' (x 3 > z+), (39) 

8 = T Vm^m?- (40) 

The coefficients b m are functions of e and r. 



We prove in section 5.3 (Theorem 16) that £ is differentiable with respect to e and r if 
these are measured in an appropriate L p norm, with p < oo, as long as e and r are bounded 
from below and above by positive numbers and there are no resonant frequencies for the 
scattering problem. The derivative is expressed in terms of the solution u a d to an adjoint 
problem in which the incident field is obtained by sending the transmitted field of u back 
toward the slab from the right. The incident and scattered fields have Bloch wavevector — k. 

V-rV?i ad + cAw ad = inf 3 \E, (41) 

M a d and T<9 n w a d are continuous on S, (42) 

M a d is — K-pseudoperiodic in (xi,x 2 ), (43) 

«ad = + <d> with < c d outgoing, (44) 



7/ mC 
"ad 



b m e- irimX3 e- i{ - m+K)x ' '. (45) 



If u + u is the solution of the scattering problem with perturbed coefficients e + e and r + f 
and ^ is the corresponding change in the transmitted energy, that is, 



£ + £ = lm / T (u + u)d n (u + u), (46) 

then the linear, leading-order, change £q in £ is given by 

£ = Im / (r Vm • Vu ad - to 2 e nw ad ) . (47) 
Jq 



It is proved in Theorem 16 that £q is bounded as a function of e and r measured in an L p 



norm (p < oo) and that the error is estimated by the square of the L p norm, 



\£ -£o\<C (||f||jy.(n) + \\e\\ L p(n)) 2 ■ (4. 
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3.2 Variation of the complex transmission coefficients 



The transmitted energy 8 is a function of the coefficients of the transmitted propagating 
harmonics b m , m G Z p . In fact, one can obtain the variational gradient of each complex 
coefficient individually. The associated adjoint problem for b m is obtained by replacing the 



incident field (45) by a single incoming harmonic, 

The linear, leading-order, change (b m )o in b m as a function of variations of e and r is 

(6m)o = [ [rVu- VuZ - u 2 e uu™) . (50) 

Compare the formulas in [8] and [10j . §4, for the case of conical diffraction by two-dimensional 
periodic structures, in which the interfaces between contrasting homogeneous dielectrics are 
varied. 



3.3 Structures with homogeneous components 



An important class of periodic structures is comprised of those that consist of homogeneous 
components. The variational gradient (47) can be formulated in terms of the material and 



geometric parameters of these components. Suppose that one period of the slab consists of 
components described by N disjoint domains Dj with material coefficients given by spatial 
constants Ej and Tj. The coefficients exterior to the components are e and r , and the 
normal vector n to dDj is directed outward. 



Variation of the values of Tj and Sj . 



If we keep the boundaries of the domains D~ fixed 



and perturb the numbers Sj and r,- by amounts e,- and r,-, then (47) becomes 



£ = Im J2\ ?j \ 



Vu ■ Vu ad - U 2 £j I ////,,! 



Im y^irj u ad d n 
3=1 V JdD > 



UU ad 



(51) 



Since the domains Dj are fixed, the estimate (48) yields 



N 



\s-s \<cJ2(\^\ 2 + \^\ 

3=1 



(52) 



and therefore (51) gives the gradient of S with respect to the numbers Ej and Tj. 

Variation of the boundaries. Let us now hold Ej and Tj fixed and let each boundary dDj 
vary in the direction of a given vector field Vj defined on dDj by allowing the points on dDj 



to flow in the direction of Vj for a distance h. Then (47) becomes 



N 

£ = Im } j j (±(r j - T )Vli • V« ad =F UJ 2 (Ej - £ )«Wad) 
.7=1 JD > 



(53) 
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in which Dj denotes the region traversed by the boundary points. The upper sign is taken 
if the normal component of Vj is directed out of Dj and the lower sign is taken if the normal 
component of Vj is directed into Dj. For small h, we obtain 

N 

£ = him / ((rj- - t )Vm • VM ad - uj 2 (ej - e )uu ad ) v y n + o(h). (54) 
4=1 J 9D 3 



From (|48|), we obtain, for sufficiently small h, 

T \v)\v r n\) . (55) 



This result implies only that 8 is differentiable with respect to h r for < r < 2/p at /i = 0; 
in particular, £ is Holder continuous with respect to uniform perturbations of the boundary. 
As discussed in the Introduction, it has been proven in two-dimensional cases that £ is in 
fact differentiable with respect to h. 



4 Eigenvalues and the scattering problem in weak form 

The weak formulation of the scattering problem places it within the framework of sesquilinear 
forms in the Hilbert space H^(Q) of /t-pseudoperiodic functions on a period Q of the scatterer. 
It allows proper treatment of guided modes, as well as existence, uniqueness, and bounds of 
solutions. The weak formulation requires the Dirichlet-to-Neumann map that characterizes 
outgoing fields. For bounded scatterers in M 3 , one may refer to Lenoir, et. ai, [15] or Colton 
and Kress [3] (§5.3); for periodic structures, our formulation is essentially the same as that 
used by Bonnet-Bendhia and Starling [1]. 



4.1 The weak formulation of the scattering problem 

By treating the Helmholtz equation in the scattering Problem[2]in the weak sense, the second 
condition on the the continuity of u and rd n u is automatically satisfied. The weak sense is 
expressed as follows: If e and r are smooth except along smooth surfaces of discontinuity, 
and if £ is a smooth vector field and h is a smooth scalar function, then a function u satisfies 
the Helmholtz equation 

V- rVu + e u 2 u = V-^ + h (56) 

at points where e and r are smooth and the condition of continuity of u and rd n u on interfaces 
between materials if and only if 

I (tVu-Vv -u 2 euv) = [ (£-Vv + hv) for all v G C^°(1R 3 ). (57) 

This weak form of the Helmholtz equation allows one to relax the regularity of e and r so 
that they are merely measurable and the regularity of £, h, u, and the distributional gradient 
of u, so that they are required only to be locally square-integrable. 
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To incorporate the pseudoperiodicity and outgoing conditions required by the scattering 
Problem [2j its weak form is posed in one period Q of the slab structure, between its bounding 
planes x 3 = z_ and x 3 = z + (see Fig. 2). The pseudoperiodicity condition is enforced by 
requiring that the solution u and the test functions v be in the pseudoperiodic Sobolev space 

Hl(n) = {ue H^tt) : u(27T,x 2 ,x 3 ) = e 2mKl u{0,x 2 ,x 3 ),u{x 1 ,2n,x 3 ) = e^ iK2 u{x u 0,x 3 )}. 

(58) 

1 

The evaluation of u on the boundary of Q is in the sense of the trace map H 1 ^) — > H^{dVt). 

The outgoing condition is enforced through the Dirichlet-to-Neumann operator for outgo- 
i _ i 

ing fields, T : H 2 (r) — ► H 2 (r). It acts on traces on V of functions in Hl(Q) and is defined 

i 

through the Fourier transform as follows. For any / £ H^(T), let f m be the Fourier coeffi- 
cients of e~ tKX /; this is a pair of numbers f m = (f~, /+), one giving the m th pseudoperiodic 
Fourier component of / on T_ and the other on r + , 

f(xi,x 2 ,z±)= f^ m+K >'. (59) 

m6Z 2 

Then T is defined by 

T-.H^^^H-kv), (ff) m = -t Vm f m . (60) 
The operator T has a nonnegative real part T r and a nonpositive imaginary part Tf. 

T = T r + iT i} (61) 

(rp~~c\ _ J W)mfm if 777> £ Z e: (PS)\ 

\ otherwise. 

Oj\~f\ — J Vmfm if 771 £ Z p , ffi^l 

\ otherwise. 

i _i 

The adjoint of T with respect to the pairing (/, g) = J r fg, for / £ Hz(r) and g £ H 2 (r) is 

r*:#i(r)->£ri(r), {T T f) m = if) m f m . (64) 

T characterizes the normal derivative of an outgoing function on T as a function of its values 
on T. If we denote the trace of u on V by w again, then 

d n u + Tu = on T for it outgoing, (65) 

whereas the adjoint T* characterizes incoming fields, 

d n u + T*u = on T for u incoming. (66) 

Using this together with the decomposition u = u mc + u sc of the solution to the scattering 
Problem |2j we obtain 



d n u + Tu = d n u mc + Tu 1 ' 



-2ir] m aZe irimX3 e iim+K)x ', x £ r_, 
-2ir] m ^ c e- inmXs e- i(m+K)x \ xeT + . 

m£Zp 
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The function (<9 n + T)w inc gives rise to an element of the space ff*(fi)* of bounded conjugate- 
linear functionals on H^(Q), which we denote by /p. We emphasize only the dependence on 
the frequency u, the parameters k, Eq, and tq being fixed. We also write for T. 

/p» = r J {{d n + T")u mc )v dA for all v G H^to), (68) 

in which evaluation of v on T is in the sense of the trace map. 

Problem 4 (Scattering of an incident wave, weak form) Find a function u G H].(Vt) 
such that 

[ (t\7u ■ Vv - u 2 euv) + r [ (T u u)v = f£{v) for all v G (69) 
Jn Jr 

The scattering problem is generalized by allowing /p to be replaced by a general element 
/ G ff*(Q)*. Problem [3] has the weak form 

Problem 5 (General scattering, weak form) Find a function u G H\(Vt) such that 

[ (rVu ■ Vv - uo 2 euv) + r I (T^u)^ = ff(v) + I (£ • + ^) /or a// u G ff*(^). 
Jn Jr Jn 

(70) 

The vector field £ and the function /i are in L 2 (f2), making the right-hand side a bounded 
conjugate-linear functional on if*(Q). The equivalence of the scattering Problems [2] and [4] 
as well as their generalizations [3] and [5] is expressed in the following theorem, whose proof is 
standard. 

Proposition 6 (Equivalence of strong and weak forms) Let e and r be bounded and 
measurable in Q, and let £ and h be in L 2 (Q). If u G Jf/ 0C (1R 3 ) satisfies the scattering 
Problem^ (resp. in which the Helmholtz equation ( 21f22 ) and the interface conditions 



(23) together are replaced by the weak condition (57), then u\q G H^(Q) satisfies Problem^ 
(resp. [5]]. Conversely, if u G H\{Vl) satisfies Problem\j\ (resp. \o^), then there exists a unique 
extension u of u to IR 3 such that u satisfies Problem (resp. wA). 

The unique extension u of the solution u of Problem |4] or [5] to all of space, mentioned in 
Proposition [6j admits the Fourier expansions ( 9fTo" ). Because of this, one can prove that u is 



bounded in any finite domain in M 3 by u and the incident field, as expressed in the following 
theorem. The theorem may be proved most elegantly using an integral representation formula 
that expresses the scattered field in a finite region to the left or right of one period of the 
structure as a bounded operator of the Cauchy data on V of the total field (see Lemma 2.1 
of [22] and the proof of Lemma 3.8 of [B] for boundedness) . We take a more direct approach 
here. 

Lemma 7 (Boundedness of field extension) Let D be a bounded domain in IR 3 and uj + 

a positive number. There exist numbers C% and C2, independent of uo as long as uo < uj + , 
such that, if u E if/ 0C (lR 3 ) is n-pseudoperiodic in x' = (x\,x%) and admits expansions j ^Jd ), 
then 

\\u\\ HH d) < CAuWrnm +c 2 J2 (KT + ICf)(i + M). (71) 
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Proof. Because u is pseudoper iodic, it is sufficient to prove the theorem for a domain Cl of 
the form 

Cl = {x G M 3 : < xi < 2vr, < x 2 < 2tt, z < x 3 < z_}, (72) 

and for a domain of an analogous form with z + < x 3 < zq. The proofs are analogous. It is 
convenient to express the form (Jol) as 



u(x',x 3 ) = ^ c e irim{xs - z - ] e l{m+K)x ' + a m e _!,,m(l3 " 2 - ) e i(m+K)l ' (x 3 < z_). (73) 



meZ p meZ 2 

Denote the first sum by u\ and the second by u 2 . 



I \u 2 \ 2 = 4vr 2 [ \a m e~ ir]m{x3 ~ z ~ ] \ 2 dx 3 = An 2 \a m \ 2 f \e- 2i ^ x *- z ^\dx 3 

< 4n 2 (z_ - z ) J2 Kl 2 + 4vr 2 ^7- ( 74 ) 

m$_Z e m£Z e 

The gradient of u 2 in Q is 

Vm 2 = a m e" ! ""' (l3 " z - ) e ,( * +m)l '(i(K + m), -i Vm ). (75) 

Similar estimates yield 





2 + 


k + m 


2 


2|r/ m 





(76) 



/ | Vu 2 | 2 < 4tt 2 (2;_ - z ) Y \am\ 2 (Vm + |k + m| 2 ) + 4vr 2 ^ |a 

For we obtain the estimate 

/ (M 2 + |V Ml | 2 ) < 4vr 2 (^ - z ) l«£T(l + Vl + (K + mf). (77) 

From the estimates ( 76|77 ) and the definition of the numbers r? m , one infers that there is a 
positive constant c such that 



c f(\u\ 2 + \Vu\ 2 )< KT(i + |m|)+ M 2 (i + H), 

^ -mez p mez 2 



(7E 



and c does not depend on uj as long as uj < uj + . The trace theorem allows us to estimate the 
coefficients a m in terms of u in Q, 

K c + a-| 2 (l + |m|) + J2 M^ 1 + \ m \) = lklr_ fei/2 (r) < M NI W)" (79) 

m&Zp ra^Zp 

From estimates ( |77|78|79[ ), we obtain 

c|M|^ ( nj < M|| M ||^ 1(C) + 2^T |< c | 2 (l + |m|). (80) 

A similar estimate is obtained for a domain analogous to Q for x 3 > z + . As a result of the 
pseudoperiodicity of u and the boundedness of D, we obtain the desired estimate. ■ 
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4.2 Eigenvalues of the scattering problem 

We will assume that the functions e(x) and t(x) are bounded from below and above in fl by 
fixed positive constants, 

< e°_ < e(x) < e° + and < r < t(x) < r° for xEVt. (81) 

The weak form of the scattering problem can be expressed in terms of the following 
sesquilinear forms in H\(Vt): 

a u (u, v) = / t Vu ■ Vv dV + r / {T^ujv dA, (82) 
Jn Jr 

a%(u, v) = / r Vu ■ Vv dV + r / (T?u)v dA, (83) 



a£(u, v)=r J (T?u)v dA, (84) 
b(u, v ) = / euv dV. (85) 



Observe that a w = + iaf. In terms of these forms, Problem [4] can be written as 

a"{u, v) - u 2 b{u, v) = fg(v) for all v G fl^(fi), (86) 
and the generalized scattering problem as 

a w («,«) -cj 2 6(u,v) = /(v) foraU«G^(fl) (/ G H*(£V)*). (87) 
We consider first the homogeneous problem 

a w (u, u) - u 2 b{u, v) = for all u G fl^(fi). 



This is a nonlinear eigenvalue problem because of the dependence of a u on uo through the 
Dirichlet-to-Neumann operator T w . 

Definition 8 A number u is said to be an eigenvalue of a one-parameter family of bounded 
sesquilinear forms c w (-, •) in a Hilbert space H if there exists a nonzero element u G H such 
that, for all v G H , c w {u,v) = 0. 

The eigenvalues of the family a u — uo 2 b are in general complex. Its real eigenvalues 
form a subset of the eigenvalues of the real part of the form, namely a" — u 2 b, as stated in 
Proposition [9] below. For an eigenfunction of the real form to be an eigenvalue of the complex 
form also, all of its propagating Fourier harmonics must vanish. This means that, as long 
as Zg is empty, a nontrivial solution of (88), for real u 2 , falls off exponentially with distance 
from the slab structure; such a field is a guided mode of the slab. If this frequency to is large 
enough so that Z p is not empty, then to is an embedded eigenvalue for the K-pseudoperiodic 
operator corresponding to the partial (in x' = (xi,^)) Floquet-Bloch decomposition of the 
Helmholtz equation in IR 3 . Typically, an embedded eigenvalue is not robust with respect 
to perturbations of k, e, or r because the condition (91) below that the coefficients of all 
propagating harmonics vanish is generically not satisfied. The existence of a guided mode 
requires special conditions, such as symmetry of e and r for k = 0. The reader is referred to 
0]j an d [21] for further discussion of non- robust guided modes in this context. 
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Proposition 9 (Characterization of real eigenvalues) If uj 2 G R, then a function 



u G H\(p) satisfies the homogeneous problem (88) if and only if it satisfies the equation 

<(w,t>) -ia%{u,v) -uj 2 b(u,v) =0 for all v G H^Q) (89) 

and if and only if it satisfies the pair 

< {u, v) - uj 2 b(u, v) = for all v G H\ (tt) , (90) 
(u\r) m = for all m G Z p . (91) 



Proof. We prove that (88) is equivalent to the pair pplj). The equival ence to equation 
(89) is proved similarly. Suppose that uj and u^O satisfy (88). The imaginary part of this 



equation with v = u, together with the expression (63) for Tj, gives 

T(t ) -Vm\(u\r)m\ 2 = / (T?u)u = 0. 



E 

mdZ v 



(92) 



Since 7] m > for all m G Z p , all propagating Fourier coefficients (u\r) m of u on V vanish. This 
in turn proves that ai(u,v) = for all v G H^(fl), so that u and u satisfy (90). Conversely, 



if (91) holds, then a^u, v) = for all v G Hl(Q) and therefore (90) is equivalent to (88). 



Proposition 10 (Real eigenvalue sequences) Given the bounds (81) on the functions 
e and r, the eigenvalues u of the family a^!(u,v) — u> 2 b(u,v) consist of the elements of a 
nondecreasing sequence of positive numbers {ujj(e, t)}^1 1 that tends to oo and their addi- 
tive inverses. The eigenvalues of the family a w (M,f) — u 2 b(u,v) consist of a subsequence 



(e,r)}^ =1 of this sequence, where N is a nonnegative integer (perhaps 0) or infinity. 



Proof. The proof follows [4J. As the family a^(u, v) —uj 2 b(u, v) depends only on uj 2 , we shall 
consider only nonnegative values of to. Let uj > be given, and let us consider the set of 
numbers A such that there exists a nonzero function u G Hl(Q) such that a^(u, -) — Xb(u, ■) = 
0. According to the min-max principle (see [23], §XIII, for example), this set consists of a 



strictly increasing sequence of positive numbers {X - (e , t)} c *L 1 defined by 



X-(e,r) 



sup inf 

-i<L 2 (c)«e (v^- 1 ) ± \{0} 



b(u, u) 



sup inf 

y^-i<L 2 (n)«e (^- 1 ) ± \{o} 
u e Hl(n) 





Vu 


2 dV- 


f r j r {T?u) u dA 




u 


2 dA 



(93) 



in which the supremum is taken over all /c-dimensional subspaces V k of L 2 (Q), for k = j — 1, 
and "_L" refers to the orthogonal complement with respect to the norm b(u, u) in L 2 (Q). One 
can prove that, for each positive integer j, A^(e, r) is a continuous and nonincreasing function 
of uj > (see the proof of Theorem 3.3 of jl]). There is therefore, for each j, exactly one 
positive number, which we denote by ujj(e,r), that satisfies 



X- A£ ' t) (s,t) 



(94) 
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The number uj is an eigenvalue of the family a^(u,v) — u 2 b(u,v) if and only if there exists 
an integer j such that u 2 = Xj(s, r). The sequence {ujj(e, t)}^i therefore consists of all the 
nonnegative eigenvalues of the family. 

The second statement in the Proposition follows from this and Proposition M m 



Because the Rayleigh quotient in the min-max principle (93) decreases with an increase 
in e and increases with an increase in r, the eigenvalues inherit the property of monotonicity 
with respect to these functions. 

Proposition 11 (Eigenvalue dependence on e and r) Let e_ ; e + , r and r + be mea- 



surable real-valued functions on f2 that satisfy the bounds (81) and the inequalities £-(x) < 



£+(x) and T-(x) < t+(x) on Q. Then, for each positive integer j, 

UJ j(£+,T-) < UJj(£^,T + ). (95) 



5 Proof of the main theorem 

Proof of the theorem on differentiability of the solution of u and the transmitted energy 8 
with respect to e and r rely on Meyers' theorem on higher regularity of solutions of elliptic 
equations. As we have discussed, in order to apply this theorem to the solution u of the scat- 
tering problem, it is necessary to be assured that u is uniformly bounded over all admissible 
functions £ and r. The precise condition we will need is one on lower and upper bounding 
functions for these material coefficients. 

Condition 12 (Non-resonance) For a given number w 6 R, the measurable real-valued 
functions £_, e + , r_ ; and r + on Q satisfy the non-resonance condition if, for each pair (e, r) 
of measurable real-valued functions on Q that satisfy 

£-(x) < e(x) < £+(x) and T-(x) < r{x) < r + (x), (96) 

for all x E Q, uj is not an eigenvalue of the family — uj 2 b. 

This condition can be arranged if we choose the upper and lower bounding functions such 
that, for some integer j, 

Uj (e_ , t+) < u j+1 (e+, r_ ) . (97) 



Then, Proposition [TT| guarantees that, for all functions e and r between these functions, 

Uj(e, r) < LUj(e^, r+) < u j+1 (e+, r_) < w i+ i(e, r), (98) 



in which case Condition 12 holds for each ui strictly between tOj{£—, r+) and cj 3 - +1 (e + , r_). 
The condition (97) can be achieved, for example, by beginning with a fixed pair of functions 
(e, r) for which Uj(e, r) < Uj + i(e, r), and varying them up and down continuously in the L°° 
norm, with respect to which each Uj(e, r) is a continuous function of e and r. In view of the 



fact that the eigenvalues of a u — uo 2 b are a subset of {u)j}j° =1 , the condition (97) is evidently 
stronger than what is necessary. In fact, as we have discussed, any ujj is typically not an 
eigenvalue of the scattering problem, for given material coefficients e and r. 
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5.1 A uniform bound for the scattered field 



The following theorem guarantees a bound on the solution of the scattering problem that is 
uniform over functions e and r bounded below and above by functions that satisfy Condi- 
tion I 



Theorem 13 (Bound on the scattered field) Lets-, s + , r ; andr + be measurable real- 



valued functions on n that satisfy the bounds (81) and the non-resonance Condition 12 for 



all oj in some positive interval [u;_, lo + ] . There exists a positive number K such that, for each 
uj G [c<j_,c<j + ] and each pair of measurable real-valued functions e and r on Q that satisfy 



e^(x) < e(x) < e+(x) and t(x) < t(x) < t+(x) 



the generalized scattering problem (70) admits a unique solution u such that 



Mliri(n) < K \\f\\m K (ny 



(99) 



(100) 



where f G H^(Q)* denotes the general functional on the right-hand side of (70). 



Proof. We first prove that the scattering problem (70) admits a unique solution for the 
parameters given in the Theorem. Rewrite ( J87[ ) as 

(101) 



[a u (u, v) + b(u, v)) - {uj 2 + l)b(u, v) = f(v). 



Since both a and b are bounded forms in H^(Q), there exist linear operators A u and C u 
from H^(Q) into itself, as well as an element / G H^(Q) defined through 



(A^u, v) = a u (u,v) + b(u,v), 
{C"u,v) = -(uj 2 + l)b(u,v), 

(/>)=/(«)■ 



In terms of these objects, equation (101) takes the form 



{A" + C u )u = f. 



(102) 
(103) 
(104) 



(105) 



The operator A w is bijective with a bounded inverse because a u (u,v) + b(u,v) is coercive 
(recall that is a positive operator): 

Re (a w (iz, u) + b(u, u)) = / r|Vn| 2 + / e\u\ 2 + r I {T?u)u 

Jn Jn Jv 

>t° I \Vu\ 2 + e°_ I \u\ 2 > min{r°,£ } \\u\\ H i {n) . (106) 
Jn Jn 

Moreover, C u is compact because of the compact embedding of H^(Q) into L 2 (Q). By the 
Fredholm alternative, (105) (equivalently, (101 )) has a unique solution if A u + C W is injective, 
that is, uj is not an eigenvalue of the family a w — uj 2 b. But this is implied by the non-resonance 
Condition [12] which we have assumed for uj. 
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We turn to establishing a bound on this solution that is uniform over all functions e 
and r and numbers u that satisfy the hypotheses of the Theorem and all / in the unit 
ball in H\{VL)* . To accomplish this, it suffices to consider arbitrary sequences e n and r n 



of measurable functions that satisfy the bounds (99), a sequence u> n of numbers satisfying 
uj- < u n < u + , and sequences u n G Hl(Q) and f n G H*(fl)* with ||« n ||.ffi(n) < 1 and f n — > 
such that 



/ (r n Vu n -Vv - ui 2 n e n u n v) + r / (T" n u n )v = f n (v) for all v G 



(107) 



and to prove that, necessarily, Hitnllin-ffi) — > 0. We may as well assume (by extracting a 
subsequence) that there exists a number uj G [u;_,u; + ] such that uj n — ► co>. We rewrite this 
equation as 



/ (r n Vm„- Vf - £„«„?;) + t / {T^u^v = g n (v) for all u G 
Jn Jr 

in which the elements g n G H\{VL)* are defined by 

<7» = /» + " o£) / e n u n v + r [ (T w - T^K*;. 

Jn Jr 

We shall prove that g n — >■ in i?*(fi)*. 



'108) 



(109) 



We first estimate the third term in (109) 
r / (T u -T**)u n v = - 



meZ 2 



(110) 



in which, for simplicity, u nm denotes the ft- Fourier coefficient of u mn \r- It is straightforward 
from the definition of r] m to demonstrate that there exists a number c such that, for all 
m El? and all n sufficiently large, 



UJ LUr, ^ Z 

\V m -Vm \ < C\U -UJ, 

This allows us to estimate, using ||w n ||.ffi(n) < 1) 



2 1 1/2 



(111) 



-*T0 Yl (Vm - Vm)UnmV m 



m£Z 2 

This proves that 



< T c\u 2 -ul\ 1/2 \\u n \ r \\L2(n)\\v\r\\L2(n) < t g M 2 c\uj 2 -u) 2 n \ l/2 \\v\\ H i^). 

(112) 

r J (T u - T Wn )u n v -»• as n -> oo 
uniformly over t> G H*(p,) with ||u|| = 1. 



(113) 



The second term of ( 109 ) is estimated by 









\ e n u n v 




Jn 



< e\\u 2 - c^Hlunll^n) ||«||L2(n) < M 2 e\\u 2 - ujl\\\v\\ H i {n) , (114) 
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demonstrating that 

(oj 2 — cu 2 ) I e n u n v — > as n — ► oo (115) 



uniformly over v G H\(Q) with ||t> || = 1. 

The results ( 113|115 ), together with the convergence f n — > prove the strong convergence 



of <7 n to zero in H\(VL) . 

We shall now demonstrate that there exists a function u G H^(Q), measurable functions 



e and r satisfying the bounds (99), and an infinite subset T of the positive integers such 



that the following convergences hold, restricted to indices in the subsequence T, 



«n 


u 


weak H^Vt), 


Un — 


u 


strong L 2 (f2), 


£-n ^ 


e 


weak* 




— > EU 


strong if _1 (f2), 


r n -> 


T 


G-convergence in f2 


9n -» 





strong 



(116) 



The first and second subsequence limits are due to the uniform bound on the functions u n 
in H l (VL), the Alaoglu Theorem, and the compact embedding of H 1 ^) into L 2 (Q). The 
third is due to the uniform bound on the functions e n in L°°(Q) and the Alaoglu Theorem. 
Because of the strong L 2 convergence of u n , we obtain, for each v G L 2 (Q), u n v — > uv in 
L x (f2) (for the subsequence T), and therefore because of the weak-* convergence of e n , 

/ e n M n w -> / euv for all v G D #o(fi), (117) 

</n </n 

from which we infer that e n u n — ^ weakly in ff _1 (fi) = Hq(Q)* , or, more precisely, 
that j(e n u n ) — j(su), where j is the natural embedding of L 2 (fl) into if _1 (f2) defined by 
j(w)(v) = f n wv for veHq(Q). Since this embedding is compact and the sequence e n u n is 
bounded in L 2 (Q), we have strong (and therefore also weak) convergence of a subsequence, 
£ n Un — > w G i/ _1 (fi), and by the uniqueness of weak limits, we obtain w = eu, which 



justifies the fourth convergence in the list (116). The last convergence follows from the 



strong convergence g n — > in H^(Q) and the inclusion Hq(Q) C Hj.(Q). The existence of 
the G-limit (or H-limit) r satisfying the bounds t_(x) < t(x) < t + (x) for x G Q follows 
from Theorem 2 of Murat /Tartar [19J and the discussion in the second paragraph of that 
work (p. 21). 

The divergence V-£ of a vector field £ G L 2 (Q) is the element of H~ 1 (Q) defined by 

(V-00) = / £ • Vu for u G #o( fi ); (H8) 
whose norm is bounded by the norm of £, 

||V-£|U-i ( o) < ||£|| L2( n). (119) 
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From equation (108) and items 4 and 6 in (116), we infer 



V-r n Vu n — ► j(eu) strongly in if 



(120) 



(the action of the integral over T in (108) is trivial on Hq(Q)). Because of the strong 
convergence of V-r n Vu n , the weak convergence of u n in if x (f2) and the G-convergence of r n , 
we may apply Theorem 1 of [19] to deduce that 

r n Vu n — rV« weakly in L 2 {VL). (121) 

Because of the weak convergence of u n in H l (Q), we have, for all v G Hl(Q) (for n G T) 



t ljT»u n )v = r J^u n (T»*v) — > r J^u{T"*v) = t J (T"u)v. (122) 



We can now take the limit of each term in (108) to obtain 



/ (tVu-Vv -u 2 euv) +t [ {T"u)v = for all t> G H^ifl). (123) 
By the uniqueness of the solution to this problem, which we proved above, we must have 



u = in Jf 1 (f2). Equation (108), with v set equal to u n , gives 



£%; 2 ll«n|||2(n) > ^ / £n|«n| 2 > T° 



whence we obtain 



- |#n(Wn)| 

T -(ll«n|lHim) - IKIIiW - \g n (u n )\, (124) 



(e° + uj 2 + T^)\\u n \\' 2 L 2 {n) > T^\\u n \\ 2 Hl[n) - \g n {u n )\- 



(125) 



From of the strong convergence u n — > it = in L 2 (Q) and the strong convergence of g n —> in 
ii*(fi)*, we deduce that H^nll^ 1 ^) — ¥ 0' as we se ^ ou ^ to do. We conclude that there exists 



a number K such that the solution u of the generalized scattering problem ( 70 ) satisfies 

IMk(n) < K\\f\\ Hm * (126) 



for all / G H*(Q) , for all functions e and r that satisfy (99), and for all uj G [ou-,u. 



5.2 Field sensitivity to L p perturbations 



This section contains the main theorem of this work, Theorem 15 , and its proof. The theorem 



makes rigorous the formal variational gradient, obtained in section 2.2[ of the solution u of 
the scattering problem as a function of the material coefficients e and r. The field u satisfies 

f (tVu-Vv- u 2 euv) + T I (T^v = r I (d n + T") u inc v for all v G (127) 
Jn Jr Jt 
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If we replace e, r, and u with e + e, r + f, and u + u and subtract from (127), we obtain 



/ (tVu- Vv-u 2 euv) +r I {T"u)v = - / (f V(« + u) ■ Vv - cj 2 e (u + u)v) . (128) 
Jn Jt Jn 

Retaining only the linear part of the right-hand side gives an equation for uq, the formal 
linearization of the perturbation of u about (e, r), 



/ (tVuo-Vv-uj 2 eu v)+tq i (T w u )v = - / (t\7u-Vv - u 2 euv) . 
Jn Jv Jn 



(129) 



The task is to prove that \\u — Uo\\ = o(\\(e, r)||) as e and f tend to zero in an LP norm. 

Recall that, if, for some vector function £ G L 2 (Q) and scalar function h 6 L 2 (Q), a 
function u G -f^o(fi) satisfies 

/ (tVu-Vv -uo 2 euv) = / (g-Vv + hv) (130) 

Jn 

for all v G i/o(fi), we say that u satisfies, in the weak sense, the partial differential equation 

V-tV?i-cAu = V-i + h in Vt. (131) 



The proof of Theorem 15 requires the following specialization of the theorem of Meyers 



on the higher integral regularity of solutions of elliptic differential equations^] 

Theorem 14 (Meyers regularity) Given a bounded domain D C M. 3 , a real-valued mea- 
surable function t(x) in D, and positive real numbers t_ ; t + , and Rq such that 

< r_ < t(x) < r+, (132) 

there exists a number Q with 2 < Q < 6 such that, for each q satisfying 

2 < q < Q < 6, (133) 

there exist constants G\ and C2 such that, given 

u G H l (D), i G L q (D), h G L 2 (D), (134) 
V -rVu = V-£ + h %nD, (135) 
B y {2R) C D, with y eD and R > R , (136) 

the following inequalities hold: 

\\Vu\\ L <,(B y (R)) < C 1 [\\u\\ L 2 {Byi2R)) + \\^\\ Lq{By(2 R)) + \\h\\L2(B v (2R))] , (137) 
\\u\\Li(B y (R)) < C 2 \\u\\m{B y {R))- (138) 



1 In Theorem 2 of |17j , we fix pi = 2, r = 2, and the dimension n = 3. The p in Meyers' theorem 
corresponds to our q here. We also enforce Q < 6, which guarantees r* > q for all g < Q, because r* = 
(r -1 — n -1 )^ 1 = 6 > Q > q. We may use equation (49) from the theorem because q > 2 > 6/5 = 2n/(n + 2) 
with n — 3. 
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The second statement (138) is a result of the compact embedding of H l (B y (R)) into 
L q (B y (R)) for q < 6 (see Theorem 7.26 of [T2] , for example). 

Theorem 15 (Field sensitivity to LP variations) Lete_, e + , r_ ; and r + fre measurable 



real-valued functions on Q that satisfy the bounds (81) and £/ie non-resonance Condition 12 
for all uj in some positive interval [uj^,uj + ]. Assume additionally that ^2 m& z p (\ a m C \ 2 + l^m°| 2 ) 
is bounded uniformly for uj G [u;_,u; + ]. Then there exist real numbers C andp>6 such that, 
for all uj G [a>_ , £*;+] and all measurable functions e, e, r, and f on Q that satisfy 



£-(x) < e(x) < £+(x) and T-(x) < t(x) < r + (x), 
£-(x) < (e + e)(x) < £+(x) and r_(x) < (r + f)(x) < t+(x), 



(139) 
(140) 



the following statement holds: 



If u G H^(Q) is the unique solution of the scattering problem guaranteed by Theorem 13 



(that is, u satisfies (127) for all v G Hj.(Q)), u + u is the unique solution of the scattering 



problem with r replaced by r + f and e replaced by e + e in (|127|), and uq satisfies the 

(141) 



approximate equation (129) for all v G H\(VL), then the linear operator 

v2 



{L P {Q)) Z -> H\Q) :: (r,e)^u 
(restricted to (e,f) admissible by ( 140[ ) ) satisfies 



and 



WuqWh 1 ^) < C (H^IIlp^) + ||e||ip(n)) 
\u - UQ\\ H i {n) < C (||f||iP(n) + ||e||ij.(n))' 



Moreover, the derivative (141) is Lipschitz continuous, that is, if(f,e) 
Uq denote the derivatives at (e^r 1 ) and (e 2 ,r 2 ), respectively, then 



(142) 
(143) 

Mo and (r, e) i— > 



o _ ^oIIh 1 ^) < C (ll rl_ ' r2 |Uf(n) + Ik — £ 2 |Up(n)) (H'rllLP(n) + ||£|Uf(fi)) , 



(144) 



so long as the functions e 1 , e 2 , r 1 , r 2 , e l + e, e 2 + e, r 1 + r , and r 2 + f satisfy the bounds 



Because the conclusion of the theorem holds for each q > p if it holds for a given p, the 
condition p > 6 could be logically be replaced, equivalently, with the condition p > 1. We 
have used the number 6 because, in the proof, p arises as a number greater than 6. 

Proof. Let e, e, r, and r satisfy the bounds in the Theorem, and let uj G [cg>_,cg> + ] be given. 
Let -B be a ball of radius R containing Q, and let B 1 and B 2 be the balls whose centers 
coincide with that of B and whose radii are 2R and 4i2 , respectively. Let Q be as provided 



in Theorem 14 for the domain B 2 and the constant Rq- Let s, q, and p be such that 

2 < s < g < Q, (145) 

2 _1 . (146) 



q' 1 +P' 1 



-l 



and s 1 + p 1 
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Let the constants C\ and C 2 in Theorem 14 be valid for both q and s in place of the q in 
the Theorem. Denote the solutions of the scattering problems in M. 3 corresponding to u and 
u + it by these same symbols. 

Theorem 13 and Lemma [7] together provide a number K\ , independent of the choice of 
e, e, r, and f, such that 



\u + u\\m(B 2 ) < K i an d \\u\\m(B 2 ) < K i- 



(147) 



Because the condition on the incident field stated in the Theorem makes /p bounded uni- 
formly over u), this number K\ is also independent of u) G [u;_,u; + ]. 

We begin by bounding Hwoll^rn) by a multiple of ||r||iP(n) + H^IUpfn)- To do this, we 
must estimate the right-hand-side of (129) in H*(fl) . Since u satisfies the scattering problem 
(127), u satisfies the differential equation 

V-tVu + u 2 eu = inM 3 . (148) 



Applying Theorem [14] to this equation yields 

||w 2 <MIl2(^) < u; 2 ||e||L P (n)||w||z,= (n) < C 2 u; 2 ||e||LP(n)||w||tfi(n) < C 2 K x u)\ \\e\\ L p^), (149) 



\tVu\\ L 2( U ) < \\f\\ L p(n) 



< Ci||r||LP(n) + \\u 2 eu\\ L 2 (Bi) ] < C x Kx{l + e 



5-^)11^11^(0). (150) 



From these estimates and Theorem 13, we infer that ito is the unique function in H^(Q) that 
satisfies equation (128) for all v G Hl(Q) and that there is a constant K 2 , independent of e, 



e, r, r, and ui such that 



|wo||fl-i(n) < K 2 \\\t\\lp(si) + \\e\\LP(n)] ■ 



(151) 



Because of this, the linear functional (e,r) 1— > u is uniformly bounded from (L P (Q)) 2 to 
ff*(0), proving the first part of the Theorem. 

An analogous argument can be applied to the system ( 128 ) for u and the corresponding 
differential equation for u + u, which appears on the right-hand side of (128), 

V-(t + t)V(m + u) + u 2 {e + e){u + u) = in R 3 . (152) 

This results in the inequality ( 151[ ), with u in place of Uq, which, together with Lemma [7j 
yields (reusing the constant K 2 ), 



m\m(B 2 ) 



< Ko II T 



\LP(Q) 



(153) 



Next, we bound the "second-order" part of the right-hand-side of (128), namely ui eu and 



V • tVu in H^(Q)* by a multiple of (||f ||.z>(n) + • For this, we apply Theorem 

to the differential equation 



V-tVu + eu u 



-V-fV(« + u) — eu (u + u) in B 2 



(154) 
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||co> 2 eM|| L 2 (n) < ^ 2 ||£||LP(n)||«||L=(n) 

< C 2 uj 2 \\e\\ L p(n)\\u\\m(n) < C 2 K 2 uj 2 + \\e\\ L p^) [\\t\\lp{0) + ; (155) 

\\rVu\\ L 2 {n) < \\T\\ LP(n) \\Vu\\ L s {n) 

< Ci||-r||iP(n) [INU 2 (Bi) +^ 2 |k^lU 2 (Bi) +uJ 2 \\e(u + u)\\ L 2 {Bl) + ||rV(« + u) \\l^b 1 )] • (156) 
The first two terms in the right-hand-side of this last estimate are in turn estimated by 

||m|U2 (Si ) + u 2 \\eu\\ L 2 {Bl) < K 2 (l + uj 2 + e° + ) [\\r\\LP(n) + PIU*>(fi)] • (157) 
The last two terms are estimated by 

ui 2 \\e(u + w)||l2 (Bi) < cj 2 P|| LP ( Bi) ||m + m|| Ls(jBi) 

< C 2 uj 2 \\e\\ LP ( Bl )\\u + < C 2 i^io;^||e||i P ( Bl ), (158) 



\\tV(u + u)\\ L s( Bi) < \\t\\lp(b 1 )\\^{u + u)\\li{b 1 ) 
< Ci\\f\\ LP{Bl) [\\u + u\\ L 2 {B2) + \\u 2 (e + e)(u + u)\\ L 2 {B2) ] < dKi(l + £ + u 2 + )\\t\\ Lp{Bi) . 

(159) 

To simplify notation in the rest of the proof, the symbol C will denote different constants. 



By putting the estimates (156, 157, 158, 159) together and using the periodicity of e and r 



to get [||r|| LP(Bl) + ||£|Up(b7)]< C [||f||ip(n) + ||e|Uf(n)], we obtain 

||^Vm|| L 2 (q) < C||f||ip(n) [||T|| iP( n) + ||e||n>(n)] • (160) 
Considering L 2 functions and their divergences as elements of if 1 (f2)*, we conclude from 



(155) and (160) that 



V-rVw + io 2 eu\\ H1{n)t < C (||r|| L p(n) + ||£|Up(q)) 2 • 



(161) 



Equations (128) and (129) give the system for u — u , 



/ (r V(u — Mo) • — uj 2 e (u — uo)v) + r / T u {u — uq)v 
Jn Jr 



(tVu-Vv -u 2 euv) for all v G H^(Q). (162) 

Jn 



This, together with Theorem 13 and equation (161), gives us the desired result 

, 2 



\u - UqWh^O) < C (||r||it.(n) + ||£||z,p(n)) , 



(163) 



in which the constant is independent of e, s, r, f, and u, subject to the conditions in the 
Theorem. 
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Finally, we prove the Lipschitz continuity of the derivative with respect to e and r in the 
L p norm. To do this, we let e and f be fixed as directions of differentiation and perturb the 
functions e and r, at which we differentiate, by functions e and f , remembering to require 
that s, t, e + e, r + r, e + e, r + f, e + e + e, and r + f + f satisfy the bounds in the hypotheses 
of the Theorem. This results in perturbations of the fields u and uq, 



(e,r) 
u 

Uq 



[S + £, T 
U + U, 
U + U . 



T 



(164) 
(165) 
(166) 



By subtracting (129) as it is written from (129) with these substitutions, we obtain 



((r + f)Vuo • - lo(e + s)u v) + r / T^uqv 
o Jr 

— — (rVn ■ Wv — uj 2 suv) 



(tVu • Vv 



UJ 2 SUqV) 



Steps analogous to those between equations (155) and (161) give us the estimates 



I V-tVm + uj 2 su\ 
I V-rVwo + u 2 eu \ 



Ih 1 ^)* 



< c 

< c 



I^IUp(n) + H^IUf(n)) (||^IUp(n) + ||^IUp(n)) , 



(167) 



(168) 
(169) 



in which the constant is independent of the functions. An application of Theorem 13 



to 



equation (167) gives a uniform bound 

H^oH/f^n) < C (||f||xp(n) + ||£||LP(n)) (||^||LP(n) + 



(170) 



This means that the linear functional (e, f) i— > {uq + Uq), which is the derivative of the total 
field with respect to the material parameters at (e+e, r+f), is bounded by C (\\t\\lp(q) + ||e||.LP(n)). 
We conclude that the derivative, defined on admissible functions e and f , is Lipschitz con- 
tinuous with respect to the L p norm over functions e and r that satisfy the hypotheses of 
the Theorem. ■ 



5.3 Transmitted energy 

We take b mc = so that there is a source field incident upon the slab only from the left. The 
energy transmitted to the right-hand side of the slab is given by 

£ = Im T ud n u = —lm / t (Tu)u. (171) 

Let u be a arbitrary perturbation of u, and 8 the corresponding perturbation of 8, 

8 + 8 = -lm / T (T(u + u))(u + ii). (172) 
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From the equations for £ and £ + £, we obtain 

£ = -Im / t ((Tu)u + {Tu)l + (Tu)u). (173) 

Denote the linear part of £ by £q . 

£ (u,u) = -Im / to((Tm)h + (T«)«) (174) 



Because of the trace theorem and the boundedness of T, we have 

£<C\\u\\ 2 H1(n) , (175) 

|£| < C\\u\\Hi(n)\\u\\m(n) , (176) 

|£-£>| <C||«||^ 1(n) . (177) 

This demonstrates that the map H^(Q) — > R :: u i— > £ is bounded and different iable and 
that the bounded linear map 

i^(fi) -»• R :: « t-> £ (178) 



defined through (174) is the derivative oi u ^ £ aX u. 

Because the adjoint problem for the transmission described below is a scattering problem 
with wavevector —k, we will need to exhibit explicitly the dependence of rj m , Z p , and T on 
the wavevector, which we have suppressed until now. From the definition (J8]) of r\ m and Z p , 
we obtain 

= Z~: = -Z;. (179) 



Using this and the relation = u K _ m for the k- Fourier coefficients (59) of a function u G 
H^(Q) restricted to V by the trace map, one can derive the relations 



T*u = T_ K u , T K u = T* K u , (180) 
and thence the equivalent expression for the differential of the transmitted energy 

£ = £ (u,u) = -Im [ T {{X K -T* K )u)u. (181) 



Essentially following [16] . we will demonstrate that the solution w a d to an adjoint scat- 
tering problem represents the transmission functional. We take as incident field a — k- 
pseudoperiodic left-traveling wave w^f, incident upon the slab from the right, that is ob- 
tained by sending the transmitted propagating harmonics of u back toward the slab. This 
is done by conjugating u and retaining only the propagating harmonics: 

u = h m e ir] ^e l{m+K)x ' {x 3 > z+), (182) 
= propagating part of u = ^ b. m e- iri ^ X3 e i{m ' K)x ' '. (183) 
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Let M a d £ be the solution to Problem [4] with Bloch wavevector — k and incident field 

u™', thus -u a d satisfies 

(rVu^Vw-cuhu^+To [ {T_ K u &A )w = rj (d n +T_ K )u™w for all w G H} K (Sl). (184) 



Observe that, since u 1 ^ has no right-traveling part, (d n + T. K )v%£ \ r — (T_ K — T_* K )u™\ v+ By 
making the identification uq = w, and using the first identity in (180) and the definition of 
M ad C together with the fact that (T_ K — T*) vanishes on the linear and exponential Fourier 
harmonics, we obtain 

/ (rVu ■ V^ad - ^suqu^) + t {T K u )u ad = t / ((T K - T*)u)m . (185) 
Jn Jr Jr + 

Using the identification w a d = v in equation (129) for the derivative of (r, e) \— > u and 
equation (181), with uq in place of u, we obtain an expression for the derivative of the 



composite operation (r, e) i— »■ z 

4 = £o(«, Wo) = Im / (f Vu ■ V« ad 
./n 



)• 



(186) 



The derivative of (r, e) i— ► w i— ► £ is Lipschitz continuous, as we now demonstrate. 
Consider the derivative at two different pairs (e,r) = (e 1,2 ; r 1 ' 2 ), 

(187) 



(e, r) I— >■ Mq h- > £(w\ z*q) at (e 1 ^ 1 ) 



(e, r) i— > Uq i— > £(w 2 , 



at 



(e 2 ,r 2 ) 



By using equation (181), we obtain 
£{u\ul) - £{u\ul) = -Im / r 4(T K -r:)(tl 1 - M 2 )^ + (T K -r:)il 2 (^-^)]. (189) 

From of the boundedness of (T_ K —T_* K ) and estimates ( 153|151|147|144 ), we obtain the estimate 

\£{u\ul) - £{u\ul)\ < C (H^-r 2 !!^) + ^ -e 2 \\ LP(n) ) (||f \\ LP(Q) + \\e\\ LHQ) ) , (190) 

which demonstrates the Lipschitz continuity. 

The results of this section are summarized in the following theorem. An analogous 
theorem can be established for the variational derivative of b m , given by (50). 

Theorem 16 (Transmission sensitivity to L p perturbations) Let the stipulations in 
Theorem 15 hold, as well as the equations (171\172 ) defining the transmitted energy £ and 
its perturbation £. Then the linear operator 

(L q (Q)) 2 — > R :: (f, e) h-> £q (191) 



defined through (186) is bounded and 



\£-£ n \ <C \\f\ 



+ ||e||LP(n) ) ■ 



(192) 



Moreover, the derivative operator (191) is Lipschitz continuous, restricted to functions (e,r) 



that satisfy the hypotheses of Theorem 15 
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